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Abstract 

Using the geometric definition of the topological charge we decom¬ 
pose the path integral of 2-dimensional U(l) lattice gauge theory into 
topological sectors. In a Monte Carlo simulation we compute the av¬ 
erage value of the action as well as the distribution of its values for 
each sector separately. These numbers are compared with analytic 
lower bounds of the action which are relevant for classical configura¬ 
tions carrying topological charge. We find that quantum fluctuations 
entirely dominate the path integral. Our results for the probability 
distribution of the Monte Carlo generated configurations among the 
topological sectors can be understood by a semi-phenomenological ar¬ 
gument. 


* Supported by Fonds zur Forderung der Wissenschaftlichen Forschung in Osterreich, 
Projects P11502-PHY and J01185-PHY. 




The last few years have seen an ever increasing interest in studying topo¬ 
logical ideas in lattice gauge theories (see ||I|] for a recent review). Many of 
these studies were motivated by computing statistical properties of instan- 
tons such as their average size and density distribution which are of interest 
for semiclassical models p (see [Q for an analysis of these ideas for the case 
of gauge group U(l) and 2 dimensions). Rather than using the lattice to 
extract information on some underlying classical conhgurations, in a series 
of articles 0 we directly investigate possible implications of the existence of 
a topological structure for the lattice path integral. 

Here we study the characteristic values for the action of quantum fluc¬ 
tuations and compare them to analytic lower bounds of the action in each 
topological sector. Classical configurations carrying topological charge typi¬ 
cally have action in the vicinity of these bounds. It is interesting to analyze 
whether after adding quantum fluctuations the structure of the classical con¬ 
figurations is still relevant for the physical picture. 

The model under consideration is U(l) lattice gauge theory in 2 dimen¬ 
sions. Due to the simplicity of this model good statistics is easy to obtain. 
The model has the further advantage of possessing a simple form of the 
topological charge; the geometric definition based on Liischer’s work |p can 
be computed easily without making use of sophisticated techniques such as 
cooling or analyzing the spectrum of the Dirac operator. These two facts 
allow to explore the physical mechanisms governing the role of topology on 
the lattice in a rather playful way. 

Before we start developing our results we hnd it useful to comment on 
the relevance of topological arguments in the continuum and on the lattice. 
In the continuum it is possible to classify classical, i.e. differentiable conhg¬ 
urations with respect to their topological charge. It is crucial to note that 
these classical conhgurations are of measure zero in the continuum path in¬ 
tegral @] . Thus in the continuum topological arguments cannot go beyond a 
semiclassical analysis. On the lattice the situation is diherent. In particular 
for the case of U(l )2 lattice gauge theory one can assign a topological charge 
to all conhgurations (except for so-called exceptional conhgurations which 
are of measure zero). For SU(N )4 it is believed, that the path integral in the 
continuum limit is dominated by smooth conhgurations which also can be 
assigned a topological charge. Thus in a certain sense, the lattice formulation 
is more powerful for studying topological ideas and the lattice path integral 
can indeed be decomposed into topological sectors. In particular the outlined 
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questions can be formulated in a meaningful way and a consistent physical 
picture can be obtained. 


We work on a 2-dimensional square lattice A with length L. Lattice sites 
are denoted as x = (xi,X 2 ) with x* = The lattice spacing is set 

equal to 1. The gauge helds are group elements U^{x) G U(l) assigned to 
the links between nearest neighbours x, x -1- /i and their action is given by 


5 = ^9^ [l - Ref/p(x) 

xGA 


( 1 ) 


The plaquette element is dehned as Up{x) = Ui{x)U 2 {x + l)Ui{x -|- 2) f/ 2 (x). 
The gauge helds obey periodic boundary conditions f/^(L-|-l, X 2 ) = U^{1, X 2 ), 
t/^(xi,L 1) = t/'^(xi, 1). 

The model was simulated with the hybrid Monte Carlo algorithm!]^ (this 
calculation is part of a another project including also fermions IQ) with 10- 
step trajectories adjusted such that the acceptance rate in the Monte Carlo 
step was 0.8 in the mean. For all lattices sizes (up to 16 x 16) we averaged 
10® measurements. The individual measurements have been separated by 10 
updates (for small values of {3) up to 500 updates (for /? = 6) in order to 
reduce correlations. 

To check our algorithm we compare our results for simple observables 
such as {S/L3) with the outcome of the analytic solution for the model with 
open boundary conditions (see e.g. |^). We hnd excellent agreement (within 
the small statistical errors), and only for the L = 4 lattice we observed small 
deviations due to the different boundary conditions. 

For the computation of the topological charge we use the geometric def¬ 
inition which is based on Liischer’s idea of associating a principal bundle 
to each lattice configuration and defining its topological charge through the 
topological charge of the bundle 0. The case of QED 2 was worked out in 
P. One obtains 

"W = ^ E Mx ), (2) 

where the plaquette angle 6p{x) is given by 6p[x) = —Mnf/p(x) and re¬ 
stricted to the principal branch 6p{x) G (—tt, vr). Note that conhgurations 
where Up{x) = — 1 for some x are so-called exceptional conhgurations and 
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Liischer’s definition does not assign a value z/[17] to them. Those conhgura- 
tions are of measure zero in the path integral. 

An interesting observable [|I0|, pT[| is the probability distribution|^ of the 
Monte Carlo generated conhgurations among the topological sectors. In 
Fig. |I] we show our results for a 16 x 16 lattice and (3 = 4 and 6. For 
both values of (3 we observe a symmetric, Gaussian-like distribution centered 
at z/ = 0. It becomes more peaked with increasing (3. 


-5 0 5 

V 

Figure 1: Probability distribution p{u) of Monte Carlo generated conhgu¬ 
rations among the topological sectors. The symbols give our results from a 
16 X 16 lattice for (3 = 4 and 6 (10^ conhgurations each), the curve is the 
semi-phenomenological distribution function (^). 

Naturally the following question arises: What is the mechanism behind 
this distribution of Monte Carlo configurations among the topological sec¬ 
tors? A hrst guess is, that there is a lower bound for the gauge held action in 
each topological sector. Such a bound would force the action to higher values 
when increasing \v\ and through the Boltzmann factor lead to a suppression 
of conhgurations in higher sectors. 

^ A physically more relevant observable is the topological susceptibility which is essen¬ 
tially the inverse of the width of the distribution in Fig. 1. For our purposes however, 
the whole distribution is more convenient since it contains more information than the 
susceptibility. 



V 
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Here we use two different lower bounds, a strict bound and a bound which 
holds only for sufficiently large (3 but is more interesting from a physical 
point of view. The hrst one resembles the famous result [O for classical 


Yang-Mills conhgurations in 4 dimensions. It also can be derived using the 
Schwartz inequality [|T^. Essentially the same arguments can be repeated 
on the lattice. We assume that the lattice gauge conhguration U is non- 
exceptional such that a topological charge (|) can be assigned. We use the 
Schwartz inequality to obtain 


x£A 


E 

xGA 


9p(x) 


l—Qidp{x) 


[1 


gi0p(a:)j 


< ^ [1-COS 0p(x)]x^ 
y€:A 


ep{y? 

1 —cos Opijj) 


Using (|I]) and (Q) one ends up with 


where 


27r' 


S[U] > —(5 D[U] > 


27r2 

1? 
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TT^ 


TT^ 


< D[U] 


2 


1_ ^ dpix)‘^ 

1 “ cos6p{x) 


< 1 . 


( 3 ) 

( 4 ) 


D[U] is a measure for how much the average plaquette angle differs from 0. 
This functional takes values between d/vr^ and 1. In the latter case which 
corresponds to /? —*• cxo, equals the bound valid for classical conhgurations 
on a continuous 2-dimensional compact manifold (replace by the volume 
of the manifold). 

Another formula also serves as a lower bound for sufficiently smooth con¬ 
hgurations and is more stringent. Using the fact that the function (1 — cos 6) 
(the contribution of one plaquette to the gauge held action) is convex for 
|6'| < n/2 one hnds 


1 

Z2 


E[1 


cos9p{x) ] 


> 1 — cos 



which implies 


5[f/] > pL 




( 5 ) 


It has to be stressed, that this bound holds only if for all x G A the condition 
\9p{x)\ < ti/2 is fulhlled. However for large enough P this is the case and (|^) 
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parameters 

u 

S{u)/L- 

Sph{u)/I3^ 

bound P) 

bound (1) 

L = 8 

0 

0.536(5) 

0.536 

0 

0 

(3 = A 

1 

0.556(8) 

0.553 

0.0078 

0.0193 


2 

0.62(4) 

0.603 

0.0312 

0.0768 

L = 16 

0 

0.543(3) 

0.543 

0 

0 

(3 = A 

1 

0.545(3) 

0.544 

0.0004 

0.0012 


2 

0.549(4) 

0.547 

0.0019 

0.0048 


3 

0.555(9) 

0.552 

0.0044 

0.0108 

L = 16 

0 

0.523(3) 

0.523 

0 

0 

(3 = Q 

1 

0.525(3) 

0.525 

0.0007 

0.0018 


2 

0.530(6) 

0.530 

0.0029 

0.0072 


3 

0.540(16) 

0.538 

0.0066 

0.0163 


Table 1: Data for the action in varions topological sectors. For several valnes 
of L,/9 and v we give the average valne of the action per plaquette S{h>)/L‘^ 
the phenomenological value Sph{h')/L?‘ from (^) and the values for the lower 
bounds (^) and (|^) (also normalized with L~‘^). 


holds. For small \v\ it is more stringent than the exact bound (compare 
Table 1). 

It is interesting to note, that the right hand side of (§) is the value of the 
lattice action (|^) for conhgurations which correspond to continuum helds with 
constant electric field ‘I'nvjV^ where V is the volume of the base manifold. 
In the path integral for the Schwinger model an a continuous torus was 
constructed essentially using those constant field configurations and adding 
quantum fluctuations. The lower bound (|) is thus of particular interest for 
analyzing the emergence of the physical picture in the continuum P, |T^ from 
a lattice formulation. 

In order to study the outlined idea that lower bounds of the action in 
each sector are responsible for the distribution in Fig. |l], we compute the 
average value S'(z/)/L^ of the action for each topological sector separately 
and compare this average with the values of the lower bounds (P) and P). 
Fig. p and Table 1 give our results. 

In the figure the symbols show the average value of the action per pla¬ 
quette and the horizontal lines give the value of the analytic result for open 
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Figure 2: Average value S{v)/L‘^ of the action per plaquette in various 
topological sectors. The symbols are the results of our Monte-Carlo simula¬ 
tion, the horizontal lines are the values from the analytic solution for open 
boundary conditions [7]. The curves connecting the symbols come from the 
semi-phenomenological formula (|^) to be discussed below. We show our re¬ 
sults for L = 8,16 and P = 2 (circles), [3 = A (squares) and (3 = 6 (diamonds). 


boundary conditions 0 which is a sum over all sectors. The curve connect¬ 
ing the symbols comes from the phenomenological formula (^ below. The 
average value of the action increases with increasing values of |i7| and thus 
the distribution of Fig. |l| can be understood through the suppression by the 
Boltzmann factor. 

What is indeed surprising is the fact, that the values of the analytic lower 
bounds (^) and (|) are far below the actual average values of the action in 
the path integral (see Table 1). For all sectors the action stays close to the 
analytic value p[ which is a sum over all sectors (horizontal lines). Only the 
center of the distribution slightly dips below this value. The range of values 
for the analytic lower bounds (and thus the action for the constant field con- 
hgurations [Q) is typically one order of magnitude smaller (see Table 1). 


To study the quantum fluctuations further we analyze the distribution of 
the values of the action in a hxed sector. It can be computed by binning the 
range of values of the action and counting the number of conhgurations in 
each bin. In Fig. ^we show histograms for the probability distribution of the 
action per plaquette in the sectors with i/ = 0,1 and 2. The vertical lines 
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show the average value S{i') of the action in each sector as well as the lower 
bound (^). 
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Figure 3: Probability distribution p(S') of the action per plaquette in the 
sectors i/ = 0,1,2. The data was taken on a 8 x 8 lattice at /? = 2. The 
dotted vertical lines indicate the value of the lower bound (5) and the average 
S{h')/L'^ for each sector. 

It is obvious, that quantum fluctuations keep a large portion of the con- 
hgurations high above the lower bounds (j^) and (j^). For all sectors the 
conhgurations have action close to the value given by the analytic solution 
[§] for open boundary conditions. Comparison with the lower bound (|^) 
which gives the characteristic action for classical, topologically non-trivial 
conhgurations shows the importance of quantum huctuations. This domi¬ 
nance of the quantum huctuations was already conjectured in a study of the 
model on the continuous torus P]. We remark that we obtain similar distri¬ 
butions for other values of L and /?. For increasing L the total action of the 
conhgurations tends to be more concentrated around the peak, and the whole 
distribution is narrower, as expected. The lower bounds become even less 
important. Increasing (3 makes the distribution less symmetric by increasing 
the weight of smaller values of the action, essentially without altering the 
support of the distribution. 

Although the analytic lower bounds do not seem to govern the distribution 
of the values of the action in the fully quantized theory, it is still possible 
to use topological arguments in a phenomenological way to understand the 
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distributions in Fig. From the definition of the topological charge (Q) it 
follows that the average increase of the plaquette angle when increasing \v\ 
by 1 is given by 2njL‘^. Assuming that the effect of going from the trivial 
sector to a nontrivial sector can be taken into account by adding this average 
value to all plaquette angles one finds {Op{x) are the plaquette angles of a 
configuration in the trivial sector) 

[l-cos(0?,(a;) + z/|^) = 

xGA ^ 


^(0) + ^ ^ sin0?,(a;) + 

^ xeA 


2 71^ 1/2 

L4 


-S{0)+pL^ +0(z/Vh®) 


In the second step we assumed that |z/| is small compared to and expanded 
the cosine in vjL?. The second order term has been re-expressed in terms of 
the total action. 

The sum over sin6'p(x) may be neglected for two reasons. On one hand 
in the sum over all configurations it vanishes due to the symmetry of the 
Wilson action in the plaquette angle. Furthermore (at least in order 6p{x)) 
it is proportional to the topological charge which is zero per ansatz. We end 
up with the following formula for the average value of the action 

Sphiiy) = ^(0) + [ - A(0) + (3L^] + 0{y^lL^) . (6) 


In Table 1 we give the values of for the corresponding f3, L and u. The 
curves in Fig. ^ were drawn using (^). The formula gives reasonable values 
for the average action over a wide range of /?, L and u with a slight tendency 
to underestimate the Monte Carlo results for smaller (3. This tendency can 
be understood by the fact that 217 jL? only gives the minimal increase of the 
average plaquette angle compatible with an increase of \v\ by 1. Quantum 
fluctuations seem to lead to slightly higher values of the average plaquette 
angle and thus of the action. When increasing (3 the fluctuations are more 
suppressed and the quality of (@) increases. 

We finally remark, that (|^) can be used to understand the shape of the 
distribution in Fig. |I| and to test the ergodicity of the updating algorithm. 
Raising (^) to the exponent and normalizing the result gives the Gaussian 
distribution 




C 




-S{<d)+(3L^ 




p(z/) = A^exp 
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with normalization (note that v is an integer) N~^ = exp(—(7/2). It 
should reproduce the probability distribution of the conhgurations among the 
topological sectors. The dotted lines in Fig. |^show the Gaussian (8), and we 
hnd good agreement with the Monte Carlo data. The quality increases with 
(3 and L where formulas (|^) and (0) are more reliable. Thus the probability 
distribution of Fig. |l| can indeed be understood using topological arguments. 
We conclude that our updating algorithm is ergodic. 

Using a different dehnition of the topological charge, Bardeen et al. JH 
derive a similar formula for p{v) from the analytic solution |^. In the limit 
/5 —> cx), where the two dehnitions of the topological charge are expected to 
coincide, their result approaches our formula (^. In fact, replacing S'(O) in 
by the analytic result for inhnite lattices, proportional to [1 — /i(/3)/Jo(/3)] 
(which, however, sums all topological sectors) one recovers the distribution 
given in [P . 


The most remarkable outcome of this study is the fact, that quantum 
fluctuations entirely dominate the path integral. At the values of L and (3 
where our simulations were performed, the conhgurations which essentially 
contribute to the path integral have values of the action which are typically 
one order of magnitude larger than the analytic lower bounds. This domi¬ 
nance of quantum huctuations was already conjectured for the model on the 
continuous torus in [Q. In the underlying construction [Q the path inte¬ 


gral is essentially decomposed into constant held conhgurations that carry 
topological charge and saturate the lower bound of the action, plus a held de¬ 
scribing the quantum huctuations. The observed importance of the quantum 
huctuations shows that trying to establish a physical picture where dominant 
classical conhgurations are used to understand the interplay of topological 
charge and vacuum expectation values is questionable. 

Using a semi-phenomenological ansatz we were able to obtain a qualita¬ 
tive, and in the /? —*• cx) limit even quantitative understanding of the average 
behaviour of the action. This formula explains the probability distribution 
of Fig. 1^ in terms of a topological argument and can be used to check the 
ergodicity of the updating algorithm. 

We believe, that although the technical problems are much more severe, a 
similar study can - in principle - also be accomplished for lattice Yang Mills 
theory in 4 dimensions. The possible outcome would be a better understand¬ 
ing of the relation between quantum huctuations and classical conhgurations 
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carrying topological information. 

We thank Philippe de Forcrand and Helmut Gausterer for interesting 
discussions. C.R. G. has also profited from remarks by Erhard Seiler and 
Gordon Semenoff. 
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